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Abstract. It's well known that there is a so-called exponential- type error bound for Gaussian 
interpolation which is the most powerful error bound hitherto. It's of the form \f(x) — s(x)\ < 
Ci(c2<2)"3" || /|| ft where / and s are the interpolated and interpolating functions respectively, c\, C2, C3 
are positive constants, d is the fill-distance which roughly speaking measures the spacing of the data 
points, and \\f\\h is the /i-norm of / where h is the Gaussian function. The error bound is suitable 
for x € R n , n > 1, and gets small rapidly as d — * 0. The drawback is that the crucial constants C2 
and C3 get worse rapidly as n increases in the sense c 2 — > 00 and C3 — ► as n — > 00. In this paper we 

raise an error bound of the form \f(x) — s(x)\ < c l 1 (c' 2 d)'^ y/dWfWh, where c' 2 and d 3 are independent 
of the dimension n. Moreover, c' 2 << C2, C3 << c 3 , and c[ is only slightly different from c\. What's 
important is that all constants , c' 2 and c 3 can be computed without slight difficulty. 

AMS classification:4lA05,4lA25,4lA30,41A63,65D10. 
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1 Introduction 

First, let h be a continuous radial function on R n which is conditionally positive definite of order 
m. Given data points (xj, fj),j = l,...,N, where X — {x\, . . . , xn} is a subset of points in R n and 
the f'jS are real or complex numbers, the so-called h spline interpolant of these data points is the 
function s defined by 

N 

s(x) =p(x) +J2 c j h ( x ~ Xj), (1) 
i=i 

where p(x) is a polynomial in F^— 1 an< l c 'j s are chosen so that 

AT 

Y J C ] q{x ] )=0 (2) 
i=i 
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for all polynomials q in PP l _i and 

N 

P( x i) + ^ Cjh(Xj - xj) = fj, i = l,...,N. (3) 

i=i 

Here PP l _ l denotes the class of those polynomials of R n of degree < m — 1. 

A famous property is that the system of equations (2) and (3) has a unique solution whenever 
X is a determining set for P^-i an d h is strictly conditionally positive definite. More details can 
be seen in [9]. Therefore in our case the interpolant s(x) is well defined. 

We clarify that X is said to be a determining set for Pj£_i if X does not lie on the zero set of 
any nontrivial polynomial in PP l _ 1 . 

In this paper the function h is defined by 

h{x) :=e-^ x \ /3 > 0. (4) 

This is the so-called Gaussian function. 

In [llj Madych and Nelson raise the famous exponential-type error bound for the scattered data 
interpolation of Gaussian function, as mentioned in the abstract. The computation of the constants 
ci, C2 and C3 can be found in [8]. This error bound is very powerful. However, too many data points 
may lead to a large condition number when solving the linear system (2) and (3). There will be 
a significant improvement of the ill-conditioning if a better error bound can be obtained so that a 
satisfactory error estimate is reached before too many data points are involved. This is what we are 
pursuing in this paper. 



1.1 Polynomials and Simplices 

In this paper we will use P" to denote the space of polynomials of degree I in n variables. It's well 

known that it has dimension dimPp = ^ U ' J . We will denote dimPp by N in this section. Let 

E C R n be compact. The interpolation theory tells us that if xi, . . . , xn £ E and do not lie on the 
zero set of any nontrivial q £ P", there exists Lagrange polynomials k, i — 1, . . . , N, of degree I 
defined by k(xj) = 5y, 1 < i,j < N, such that for any / G C(E), (Uif)(x) :— Yli=i f{ x i)U{x) is 
its interpolating polynomial. It's easily seen that H(p) = p for all p £ Pp and hence the mapping 
n; : C(E) — > Pp is a projection. Let 

JV 

HUH :=m a xJ2Hx)\. 

i=l 

A famous property says that for any p £ P", 



IblU < HHll max : \p(xi)\ 

l<z<N 

The compact set E discussed in this paper will be mainly an n-dimensional simplex T n whose 
definition can be found in [2]. 

It will be convenient to adopt barycentric coordinates when discussing points in a simplex. 
Suppose Ui,..., v n+ i are the vertices of T n . Then any x £ T n can be written as a convex combination 
of the vertices: 

n+i 

X = 2J CiVi 
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where ~Y^=i Cj = 1 and Cj > for all «. The barycentric coordinate of a; is then (c\, . . . , c„+i). Let's 
define "equally spaced" points of degree I to be those points whose barycentric coordinates are 
of the form 

n+l 

(ki/l, ki/l, . . . , k n+ \/l), ki nonnegative integers with fc, = I. 

i=l 

Obviously the number of such points in T n is exactly N = dimP" . Moreover, by [T], we know that 
equally spaced points form a determiming set for P™. 
The following lemma cited from pQ will be needed. 

Lemma 1.1 For the above equally spaced points ||IT;| < ( ^ . ^ ). Moreover, as n — > oo, ||IL| — > 



2/- 1 



Now, we are going to prove a lemma which plays a crucial role in our construction of the error 
bound. 

Lemma 1.2 Let Q C P™ be an n simplex in R n and Y be the set of equally spaced points of degree 
I in Q. Then, for any point x in Q, there is a measure a supported on Y such that 

p(y) d < j (y) = P( x ) 

for all p in P™, and 

21-1' 



dW\(y) < 



I 



Proof. Let Y = {yi, . . . ,i/n} be the set of equally spaced points of degree I in Q. Denote Pp 
by V. For any x £ Q, let 8 X be the point-evaluation functional. Define T : V — > T(V) C R N 
by T(v) = (S yi (v)) Vi&Y - Then T is injective. Define ijj on T(V) by $(w) = 5 x (T- 1 w). By the 
Hahn-Banach theorem, ip has a norm-preserving extension ip ex t to P . By the Riesz representation 
theorem, each linear functional on R N can be represented by the inner product with a fixed vector. 
Thus, there exists z S R N with 

N 
J'=l 

and ||2||(ijjv)* = HV'extll- If we adopt the Zoo-norm on R N , the dual norm will be the ii-norm. Thus 
\\z\\ (m)r =\\z\\ 1 = U ext \\ = U\\=\\6 x T- 1 \\. 

Now, for any p g V, by setting w — T(p), we have 

N N 

$x(p) = 5 x (T~ 1 w) = r/}(w) = ipext(w) = ^ z i w i = X! z $vAv)- 

This gives 

N 

p( x ) = ^2 z jp(yj) ( 5 ) 

3=1 

where | Zl | + • • • + \z N \ = U^T-H. 
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Note that 



\6 X T || = sup 



w e T(V) ^ RN 



|p(*)l 



< sup 



< sup 



pey max^i,...,^ \p(yj)\ 
Ui\\ max,-=i,...,jv 



pey maxj=i,...,jv|p(y i )l 

n,|| 

2/ - 1 



2/ — 1 \ 

Therefore |zi| + • ■ ■ + \zjsr\ < ( ^ J and our lemma follows immediately by (5). 



1.2 Radial Functions and Borel Measures 

Before moving on to our main result, some background for interpolation is necessary. First, the 
space of complex- valued functions on R n that are compactly supported and infinitely differentiable 
is denoted by T>. The Fourier transform of a function in I? is 

i<x£> 



0(0 = J e- l<x ^4){x)dx. 

Then a crucial lemma introduced in [5] but modified by Madych and Nelson in [lOj says that for any 
continuous conditionally positive definite function h of order m, the Fourier transform of h uniquely 
determines a positive Borel measure /i on R n ~ {0} and constants a 7 , |-y| = 2m as follows: For all 



h(x)ip(x)dx — 



[ h\<2,n J 



E (6) 

H<2m 



where for every choice of complex numbers c a , \ a\ = m, 

E E a a+f3C a Cf3 > 0. 
|a|=m |/3|— m 
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Here \ is a function in T> such that 1 — x(Q nas a zer0 of order 2m + 1 at £ = 0; both of the integrals 

Jo<iei<i 7|ei>i 

are finite. The choice of x affects the value of the coefficients a 7 for I7I < 2m. 



2 Main Result 

Before showing our main result, we need some lemmas. First, recall the famous formula of Stirling. 
Stirling's Formula: n! ~ v / 27m(f )™. 

The approximation is very reliable even for small n. For example, when n = 10, the relative 
error is only 0.83% . The larger n is, the better the approximation is. For further details, we refer 
the reader to [3] and [4]. 

1 

Lemma 2.1 Let pi = \ and p 2 = ^f- ~ ^p. Then 

V2^p k k k <k\< V2^p%k k 

for all positive integer k. 
Proof. Note that 



can be expressed by 



Now, 



1 y/2 y/3 VI y/5 



e e e e e 



1 23 3e 4s 5io 3e 
sup <-,—,—,—, , ...> = — 

e e e e e e 



implies that ^ < p\ for all k. Thus k\ ~ v^^l ■ k k < \[zltp\ ■ k k . 
The remaining part \/ 7 2/np\k k < k\ follows by observing that 

2^{-) k k k < V2^(-) k ■ Vk, ■ k k ~ hi. 



Lemma 2.2 Let p= ^. Then k\ < V^/fc^ 1 for allk>\. 
Proof. First, 

k\ ~ V2^(-) k ■ Vk ■ k k = V2^ ■ % ■ k k -\ 
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Note that j^|:fc = l,2,3,...j can be expressed by 



1 2l 3^ 4l 



Our lemma follows by noting that 



1 2« 3* 4* V3 
sup < -, — , — , — ,.••>= • 

e e e e j e 



Lemma 2.3 Let h(x) = e ^l x l ,/3 > 0, 6e i/ie Gaussian function in R n , and fj, be the 
measure defined in (6). For any positive even integer I, 

[ \Z\ l dfi{Z) < -n-a n - 2^ • ■ /3* -(l + n- 1)^ ■ (2 + -) 

for odd n, and 

f , ,1 , . n + l l + n+3 l + n-2 I , i + n-4 

/ |£|'d/i(0 < 7T— ■ n • a n ■ 2-^~ ■ p~^~ -(3* - (l + n- 2)~^~ 

for even n, where p = & and a n is the volume of the unit ball in R n . 
Proof. 

JR^ 



R e TF 



= - n-an ■ I '- 2 dr 

P Jo ejp 

= (-)T.„.a n -(2v^) ,+n - / -^dr 



tt* • n • «„ • 2™-^ • (3l ■ (l±2fl)(l±^ - 1) ... (i) /~ 
i/ n odd, 

?rt . n . an . 2 n - 1+l ■ (3i ( l -±^)\ if n is even. 



Let 6 = pi+^+2-|. Thcn 6 | < ^pfc . Thcn for odd n ^ b = xhus 

, l+n-l ,1 + Tl — 1 , i + n-3 , -i-ra+3 i+n-1 . , i+n-3 

6! < V2np~^~ • ( )- 5— = V2^-2 — 5 — . p— 2— ■ (Z + n- l)^ - . 



; , H+1 In ' + " + 2 „I l + n-l , i + r>-3 



and 

/ KI'dMO < ^ ■ n ■ a„ ■ (2 + -) ■ 2^^ • /3* • -(l + n-l) 

JR n e 

by noting that J °° -^dv <2 + \. For even n, b = Thus 

6! < V^p^ - • ( = V2tt-2 — 5 — . . (l + n-2)^—., 
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and 



If I'dMO < ^ -n-atn- 2 i± ^ ■ 0* ■ p^^ • (Z + n - 2)^ 



Our interpolation is based on a function space called native space, denoted by Ch, m - If 
X> m = ^<j) G V : J x a 4>{x)dx = for all \a\ < mj 
then Ch.m is the class of those continuous functions / which satisfy 



< c(f) h(x- y)(j>(x)^(y)dxdy (7) 



f(x)<f>(x)dx 

for some constant c(f) and all <f> in 2? m . If / S C/i, m , let ||/||/[ denotes the smallest constant c(/) for 
which (7) is true. It can be shown that \\f\\h is a semi-norm and Ch, m is a semi-Hilbert space; in the 
case m — it is a norm and a Hilbert space respectively. In this paper both the interpolating and 
interpolated functions belong to the native space. 

The function space Ch.m is introduced by Madych and Nelson in [9] and [10]. Later Luh makes 
a lucid characterization in [6] and [7]. Although there is an equivalent expression for Ch, m which is 
easier to understand, we still adopt Madych and Nelson's definition to show the author's respect for 
them. The main result of this paper is the following theorem. 

Theorem 2.4 Let h(x) := e - ^' be the Gaussian function in R n . For any positive number bo, 
there are positive constants 8q,c\,c%, and C3 independent of n, for which the following is true: If 
f £ Ch.m , the native space induced by h, and s is the h spline that interpolates f on a subset X of 
R n , then 

\f(x)-s(x)\<c 1 VS{c 2 5)^ -\\f\\ h (8) 

for all x in a subset fl of R n , and < 5 < 5o, where f2 satisfies the property that for any x in £1 
and any number < r < bo, there is an n simplex Q with diameter diamQ = r,x £ Q C fl, such 
that for any integer I with ^ < I < 2|a, there is on Q an equally spaced set of centers from X of 
degree I — I. (In fact, the set X can be chosen to consist of these equally spaced centers in Q only.) 
Here \\f\\h is the h-norm of f in the native space. The numbers 5o,c\,C2, and C3 are given by So '■= 

7== • — }= for odd n, 




'16-n 

1 1 



lin l & o, „4 33 2 7 a \ where p 3 = 123 . ^fifj l I v^F where 



A" is defined by A" 



^2 + i • n^r- ■ (na n )^ ■ 2% ■ p~^~ for odd n, 

tt^~ ■ (iia„)5 • 2~^~ ■ p~^~ for even n, with p = — as in Lemma2.3 
,where the number a n denotes the volume of the unit ball in R n . 

In particular, if the point x in 51 is fixed, the only requirement for is the existence of an n 
simplex Q, with diamQ = r, 1 £ Q C O, satisfying the afore-mentioned property of equally spaced 
centers. 

Proof. For any b > 0, let do ■— min|& , pi.&^-b* } wnere P3 — 12* ■ \fe~P- For any < 5 < So, 
there exists an integer I such that 1 < 4-1 < 2 since < 4- < 1. Such I satisfies ^ < I < an( j 

& — bo — bo — o — — 5 

b f<ll< bo- 
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For any x € 0, let Q be an n simplex containing x such that Q C f2 and has diameter diamQ = If. 
Then Theorem4.2 of [TU] implies that 

|/0c)-s(a0l<Q||/IU / Itf-aM^GO (9) 
whenever I > 0, where er is any measure supported on X such that 

/ p{y)da(y)=p{x) (10) 

for all polynomials p in P/?^. Here 



ii" 



whenever Z > 0. Be careful. We temporarily use c\, 02,03 to denote numbers which are totally 
different from the c\, 02, C3 mentioned in the theorem. 

Let Y be the set of equally spaced centers from X of degree I — 1 on Q. By Lemmal.2, there is 
a measure a supported on Y such that (10) is satisfied and 

21 - 3 
I - 1 



dW\(y) < 

The crux of our error estimate is to bound 



I = ci[ \y-x\ l d\a\{y). 

JR™ 



Now Lemma2.3 applies. For odd n, 



* = U I \^\ 2l M0 



R" 



1/2 



1 f rt + l 2i + n + 2 2l + n-l , 2l + n-3 1 

< 77 <{ 7T— • n ■ a n ■ 2~^~ • p^— . . (21 +n- 1)^— • (2 + -) 



1/2 



2l + n~3 I 1 

• A/2 + - 



< i jvr^ • (na,,) 1 / 2 ■ 2^ • p 2 ^ . (^2^)' • (21 + n - 1) 

< /ttF (nan) 1 ' 2 ■ 2^ • ^(^)'(2I + „ _ 1)^^(2 + I) J / • p' ■ 
6y Lemma2.1 

= tt^ • (na„) 1/2 ■ 2* • p 2 ^ • p< • (i)(2i + n - l) 2 ^ ■ (2 + i)V= 

{' e 



where p$ 

Pi 



Pi 



A" • pfj • i(2Z + n - l) 2£± ^w/iere A' 
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Note that the diameter of Q is y . This gives for I > n — 3 



< 



< 



21-3 



I < A"p 3 r l (2l + n- l)—^(diamQy ^ y_ ^ 

if l>n-3 

A \ v — J. j 



< A"p l 3 l- l (3l)^(-l) 



A"( P3 3^°-i 3 /y 



I - 1 

2/ — 3 
l-l 

21-3 
I - 1 



if I >n — 3 



A" (p 3 3 3/4 • J • <5(^) 3/4 j -^=—^A l ■ \ by Stirling's formula and 
I 2 o J y 7T v/ — 1 4 



5o <?< 26o 



i 1 



< A 



16^ VI - ! 

W /iere p 3 ■ 3 3 / 4 • 2 • 2 3 / 4 • 6 3/4 • S 1 ^ < 1 iff S 1 ^ < 



p 3 3 3 / 4 • 2 7 /4 . b * 



3/4 



which 



is guaranteed by 5 < So := min | feoj —4 — 53 — ^7 — 73 \ 



< A" {pj ■ 3 3 ■ 2 7 ■ bl ■ 8}^ — 1 



< A"- 



a67r 



/T67f vT^T 

{G(5}* • -^=L= where G = /5 4 3 3 2 7 6 3 and 3 = ^ 



i{G«M(i/!»0) 

/l07T V< 

-^=A/-r-^{G(5} f since ^ < Z < ^ */ I > n-3 



For even n, if I > n — 4 
I 



q / |y-z|'d|(7|(y) 



< A"p 3 ■r l {2l + n-2) 21± t- i (dimQ) 1 



21-3 
I - 1 



where A" — tt 4 ■ (na n )^ ■ 2 4 • p 4 and p 3 = 



pi 



= A"p l 3 -r\2l + n-2Y 



2i+^-4 [5l\ [ 21-3 



¥ r 1 \* ' 2/ - 3 



< A"^.(30*-<5'-(2) 



/- 1 



/ - 1 

since I > n — 4 
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< 



< 



< 



3,3 1„, 1 



< A"( P3 3W-5) 



1 



4 by Stirling' s formula 



1 



A" 



< A" 



J_ 1 (p3.2-3i«iy (i/I»0) 

l07T V< 



1 



1 

/ T67r VI 

1 



26q 
<5 



/ 16tt VZ 
1 



p 3 • 2 • 3* • 5 
L{ P3 .2i.3*-&|.*i} 

If 7 3 2 , 1 

-=(p3-2--3-^ .51} 



sznce — < I < — — 




167T Vl 

p 3 ■ 2* • 3* • &f -Si<liffS< 
1 



where 



p\- 27. 33. 6g 



, which is guaranteed by 



8<8 a 
A"- 
A"- 
A"- 



4 27.3 3 -&3' 60 



1 1 

/ 167r VI 
1 1 



{^•2 7 -33.63.£} 



'<0 
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_ {G<5}* where G = p\ ■ 2 7 • 3 3 • b\ and g - 

I6n VZ 4 

bo . , . 26 
since — < i < — — . 
" 



jL= ■ ■ VS {G8}% 

167T V^O 



Our error bound (8) then follows immediately by letting p 3 = \fej3- 12 4 . jj 
In the proof of Theorem2.4 although we require that Z > n — 3 for n odd and Z > n— 4 for ro even, 
besides Z >> 0, it usually causes no trouble because 8-^0 implies I — ► 00 and usually 8 is very small. 
In order to make the theorem easier to understand, we don't put these requirements into the theorem. 



Remark. In theorem2.4 we avoid using the well-known term "fill-distance" for scattered data 
approximation because in our approach the data points are not purely scattered. However the num- 
ber 8 is in spirit equivalent to the fill-distance in the sense that 5 — > iff the fill-distance d(Q, Y) — > 
where Y is the equally spaced centers from X of degree I — 1 on Q. Note that d(Q, Y) — ► iff I — ► 00, 
and I — > 00 iff 8 — ► 0. 

Furthermore, although the function space Ch, m is defined in a complicated way, there is a simple 
expression for it. This can be seen in Wendland's book [12] where Madych and Nelson's native space 
9 and Wu and Schaback's native space [13] are unified in an elegant way. 



3 Comparison 

The exponential-type error bound for Gaussian interpolation raised by Madych and Nelson in [11] 
is of the form 

\f(x) - s(x)\ < 0i(02«)*||/|| hl as 8 -» (11) 
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where a\ is about the same as c\ of (8), 03 = jp- and a 2 = (3 3 / 4 • e • y/2p/3 ■ \fri ■ e 2 ™ 7 ™) 4 • b\ • j n 
where p = ^ , 60 is the side length of a cube and j n is defined recursively by 

71 = 2, 7„ = 2n(l + 7„-i) if n> 1. 

These can be seen in [IT] and [8]. Note that 7„ — > 00 rapidly as n — > 00. The first few examples are 

7i = 2, 72 = 12, 73 = 78, 74 = 632 and 75 = 6330. 

This means that for high dimensions, a 2 will become extremely large and 0,3 extremely small, making 
the error bound meaningless. 

Our new approach avoids this drawback. Although we require that the centers be equally spaced 
in the simplex, it causes bo troubles at all both theoretically and practically. 
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